Typeset with jpsj2.cls <ver.l.2> 



Letter 



O 

o 

(N 



(N 

O 
O 

03 



i 

C 

o 
o 



(N 
> 
O 

(N 

in 
o 

-I— > 
o3 



i 

C 

O 

o 



X 



Electron-Phonon Mechanism for Superconductivity in Na 35 CoO 2 : 
Valence-Band Suhl-Kondo Effect Driven by Shear Phonons 

Keiji Yada and Hiroshi Kontani 

Department of Physics, Nagoya University, Furo-cho, Nagoya 464-8602, Japan 

To study the possible mechanism of superconductivity in Nao.3sCo02, we examine the inter- 
action between all the relevant optical phonons (breathing and shear phonons) and ti g (ai g + e' g ) 
electrons of Co ions, and study the transition temperature for s-wave superconductivity. The 
obtained T c is very low when the e' g valence bands are far below the Fermi level. However, T c 
is strongly enhanced when the top of the e g valence bands is close to the Fermi level (e.g., 
— 50 meV), thanks to the interband hopping of Cooper pairs caused by shear phonons. This 
"valence-band Suhl-Kondo mechanism" due to shear phonons is important in understanding 
the superconductivity in Nao.3sCo02. By the same mechanism, the kink structure of the band 
dispersion observed by angle-resolved photoemission spectroscopy (ARPES) measurements, 
which indicates the strong mass enhancement (m*/m ~ 3) due to optical phonons, is also 
explained. 
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Despite enormous theoretical and experimental effort, 
the origin of superconductivity (SC) in Nao.35Co02 • 
2/H2O (T c = 4.5 K) is still controversial. Several 
NMR/NQR measurements above T c have revealed the 
existence of prominent magnetic fluctuations with finite 
momenta, 1-3 which might remind someone of unconven- 
tional SC due to Coulomb interactions. Recently, the 
anisotropy of Knight shift had been measured in aligned 
crystals 4 or in oriented powders. 5,6 Below T c , a sizable 
change in Knight shift is observed for both H _L c and 
H || c in aligned crystals, which indicates a singlet SC. 4 
Note that no shift is observed for H || c in oriented pow- 
ders. 6 

Significant information on the Fermi surface has 
been obtained by recent angle-resolved photoemission 
spectroscopy (ARPES) measurements: They show that 
Nan.35Co02(-?/H20) possesses a single hole like Fermi 
surface composed of a\ g orbitals in Co ions, whereas 
the other two bands which originate from e g i orbitals 
are completely below the Fermi level. 7-9 They report 
that the top of e g <-like valence bands is located about 
20 ~ 100 meV below the Fermi level, irrespective of the 
fact that a LDA band calculation predicts the existence 
of Fermi surfaces of e g i bands, which are composed of 
six small hole pockets just inside K points. 10 Theoreti- 
cally, the absence of small hole pockets would be unfa- 
vorable for unconventional SC. 11 It is noteworthy that 
the "weak pseudogap behavior" observed in the uniform 
susceptibility below room temperature in Nao.3sCo02 1 
is well reproduced by the fluctuation-exchange (FLEX) 
approximation only when small hole pockets are absent, 
as observed in ARPES measurements. 12 

Besides the on-site Coulomb interaction, the exis- 
tence of a strong electron-phonon interaction (EPI) 
in Na^Co02 is indicated by various optical measure- 
ments. 13 ' 14 In addition, the quasiparticle dispersion for 
the aig band observed by ARPES measurement has a 
prominent "kink" structure at around 70 meV below the 



Fermi level, which is expected to originate from phonons 
because the energy of optical phonons is also ~ 70 meV. 
The estimated mass-enhancement factor due to phonons 
is about three. On the other hand, the bandwidth of t 2g 
bands observed by ARPES measurement is around 0.8 
eV, which is about half the bandwidth given in the LDA 
band calculation. 10 Thus, the mass-enhancement factor 
due to Coulomb interaction is about two, 15 so the to- 
tal mass-enhancement factor becomes 2x3 = 6. This 
fact would indicate that the EPI will be dominant in 
Na^CoC^. 




Fig. 1. Displacement of O ions by E\ g (shear) mode and A\ g 
(breathing) mode. 



In the present work, we analyze the EPI on the basis 
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of the d-p model for Nao.3sCo02 given in rcf. 12, and 
find that both shear and breathing optical phonons are 
strongly coupled with t 2g electrons of Co ions. If only 
the breathing phonon is taken into account, the T c for 
s-wave SC obtained by the Eliashbcrg equation is very 
low. On the other hand, a much higher T c is obtained 
when we consider both shear and breathing phonons, 
thanks to the interband hopping of Cooper pairs be- 
tween the aig conduction band and the e' g valence bands. 
This "valence-band Suhl-Kondo (SK) mechanism" due to 
shear phonons also gives a large mass enhancement. 

The structure of the C0O2 layer is a triangular net- 
work of edge-sharing oxygen octahedra. The trigonal de- 
formation of the octahedra splits the t2 9 orbitals into 
the a\g orbital and twofold degenerate e' g orbitals. Here, 
we take the basis of t 2g electrons as |cn. g ) = (\d xy ) + 
\dy X ) + \d zx ))/VS, {eg 1 } = (\d yz ) - \d zx ))/V2, and |e^ 2 ) = 
(2\d xy ) — \d yz ) — \d zx ))/V6. We introduce V{ to express 
the change of the crystalline electric splitting between ai g 
and e' g orbitals from the value given by the LDA band 
calculation. 12 In the point charge model, V[ decreases 
and e^-orbital level is raised when the thickness of C0O2 
layer is reduced by trigonal distortion. 

First, we focus on the zone center phonon modes. 
NaCo02 has 7 kinds of irreducible representations for 
optical phonons per C0O2 layer. 16 These phonons are 
separated into a soft group (3 modes) with frequencies 
below about 400 cm -1 and a hard group (4 modes) with 
frequencies between 450 and 650 cm -1 . The latter cor- 
responds to the C0O2 octahedron's oscillating modes, 
which can give a strong EPI. In these 4 modes, we find 
that only Ei g and Ai g modes strongly couple with tig 
electrons, whereas the other two ungerade modes have 
no coupling with respect to the linear displacement of 
O ions. In Ei g and Ai g modes, O ions oscillate in the 
direction parallel and perpendicular to the C0O2 layer, 
respectively, as shown in Fig. 1. Hereafter, we call them 
" shear- mode phonons" and " breathing- mode phonons" , 
respectively. Note that the shear mode belongs to a two- 
dimensional representation. 

Hereafter, we calculate the strength of the EPI via 
the frozen phonon method. For this purpose, we ignore 
the effect of trigonal distortion for simplicity of calcula- 
tion. The displacement due to an m-mode phonon (m = 

BR, SHI, SH2) is expressed as u m = ^/ 2 AfL (^™ + ^m); 
where M is the mass of an O ion, u m is the frequency 
of oscillation and b$ is the annihilation (creation) op- 
erator for the phonons. The Hamiltonian which denotes 
EPI for a multi orbital model is written as follows: 

Hepi = ^w £ E 4 +q ^ m cka(& ro , q +&l,_ q ), (1) 

* m k,q,cr 

where c£> = (c£L 19 ,,t c^ 1(7 is the annihilation 

(creation) operator's column (row) vector for electrons, 
and V m is the coupling for m-mode phonons. Hereafter, 
we drop the q-dependence of V m and use the matrices 
for zone center phonon for simplicity of calculation. 
By choosing an appropriate coordinate axis for shear 



phonons, the matrices have the following form: 17 
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First, we consider the EPI originating from the change 
of the Coulomb potential for t 2g electrons due to the 
displacement of O ions, by considering O ions as point 

charges. The obtained result is af = — C(l — |^r), a 2 = 
-C(l + f b? = Cf and b§ = jfeb? on the order 

of a -4 . 17 Here, a is the lattice constant between a Co 
site and O site, is the ionic radius of Co 3+ , and C — 

^v^' 2 ^ vSr = °- 043 A if we p ut " = 60 

meV. By using the values a w 1.9 A and w 0.61 A , 
C is estimated to be 0.87 eV. 

Next, we consider the EPI originating from the change 

of the transfer integrals between Co and O. According to 

4 

pd-rr^^-m ' ^O-Co, 
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Harrison's law (t oc a ), St 
where eo-co = (ro — r Co)/«- The 2p orbitals of oxygen 
are filled with electrons. When we consider the virtual 
process in which a hole of t 2g orbitals transfers to 2p 

t 2 

orbitals and turns back, its energy shifts by — -P^ 1 per O 
ion, where A p d is the charge transfer energy. Thus, the 
variation in electron level accompanied by the change in 
Up to the fourth-order processes, 

-2T, and bj = T,b^ = T/y/2, 
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we obtain a\ - 
where T is given by 
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Here, we considered only nearest-neighbor hoppings be- 
tween Co— O and O—O. Using the values of Slater-Koster 
parameters given in rcf. 12, T is estimated to be 0.158 
cV. 

Thus, we obtain bi = 6f + b\ = 0.21 eV, and 
b 2 = b 2 + b 2 = 0.12 eV for shear phonons. As regards 
breathing phonons, we need to consider the chemical po- 
tential shift S[i caused by the Ai g frozen phonon since 
Tr{V BR } is finite as shown eq. (2). Because of the conser- 
vation of electron number, pd(0)(5ed + Sn) + p p (0)dfj, = 0, 
where Pd(0) and p p (0) are the density of states (DOS) 
of G?-electrons and p-electrons at the Fermi level. Thus, 
the effective shift in the 3d level is 5id = oed + Sp, = 
Pp ^ — -5ed- According to the band calculation 10 or 

MO) 



0.2. In 



pd(o)+ Pp (or 

in the present tight-binding model, Pd (o)+ P (0) 
conclusion, ai = 0.2(af + af ) <~ —0.2 eV w a 2 . Note 
that such a screening effect due to 5p is absent for shear 
phonons, so the estimated values of bi and b 2 are reliable. 
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In what follows, we analyze the strong-coupling Eliash- 
berg equation numerically. According to first-principle 
linear response calculations, 16 the frequencies of shear 
and breathing phonons are about 500 cm -1 and 600 
cm -1 , respectively, and they are almost dispersionless. 
To reduce the number of model parameters to simplify 
the analysis, we take u>d = 550 cm -1 ~ 68.2 meV for 
both phonons. In the present numerical calculation, wc 
put 



BR 



oi = 2a 2 = -a"", 61 = 4\/2&2 = a SM . (4) 

We assume a BR — a SH — 0.2, considering that the es- 
timated values of a\ and b\ are about 0.2. Both a\ and 
61 are important parameters for the SC because they are 
related to the a\ g orbital, which makes the large Fermi 
surface around the T point. On the other hand, 0,2 and 6 2 
are less important when small hole pockets composed of 
e' g orbitals are absent. We use a one- loop approximation 
to calculate the self-energy Y,(iu) n ). 
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G(k,ie n ) = r(G(°)(k,ie n )) * - t(is n ) 



(7) 



where w n = 2nnT, and e n = (2ra+l)7rT. These equations 
are solved self-consistently. 

Figure 2 shows the V{ dependence of the inverse of 
the renormalization factor z^ 1 = 1 — ^S^u) |^^ at 
T = 0.005 eV, where Y,g is the normal self-energy for the 
£ orbital. Here, zj 1 is the mass-enhancement factor for 
£-orbital because the k-dependence of the self-energy is 
absent in the present model. We have verified that z~[ Y 
is almost invariant against temperature in the present 
calculation. At (a SH ,a BR ) = (0.2,0.2), zj 1 for the a lg 
orbital decreases as V( increases. Note that the e' g band 
is completely below the Fermi level for V/ > 0.06 eV. 
To investigate the role of the two kinds of phonons, wc 
perform the calculation for (a SH ,a BR ) = (0,0.2) and 
{a SH ,a BR ) = (0.2,0) for comparison. At V{ = 0.06 eV 
where the top of the e' band is just at the Fermi level, 



about 65% of the mass enhancement for the 



tig 



band 



is brought about by shear phonons, which mediate the 
interband process, and this ratio gradually decreases as 
V t ' increases. 

Figure 3 is the renormalized band structure 
for the T-K direction in the Brillouin zone obtained 
by det|G« _1 (££) + G a ~\e^)\ = 0. At \E^\<u D , 
the electronic band is renormalized significantly and 
a kink structure appears in the a\ g band. This kink 
structure agrees well with the observation by ARPES 
measurements. 7 The e' g band is also renormalized by 
phonons, and its upper side has a flattened structure 
for \E^\ <lu d . The inset shows the top of the e' g band 
measured from the Fermi level, A, which is different from 
the original value without interaction. The small pock- 
ets disappear at V{ ~ 0.1 eV without interaction, 12 but 
they disappear at V{ ~ 0.06 eV when the EPI is taken 
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Fig. 2. Obtained 2. 1 for a\ g (left) and e' (right) orbital. 



into account. This boundary value of V{ does not change 
if only breathing phonons are taken into account, which 
means that the shear phonons depress the e' g band to a 
lower energy due to the off-diagonal elements of the self- 
energy. Although there is a linear relationship between A 
and V{, its gradient changes at V{ ~ 0.06 eV, where the 
top of the e' g band just crosses the Fermi level. We note 
that the second kink structure around 2ud is a famous 
artifact of dispersionless Einstein phonons, so it will dis- 
appear when more realistic optical phonons with small 
dispersions are assumed. 




Fig. 3. Renormalized band structure (solid line) for T-K direc- 
tion and original band structure £k (dashed line) of tight-binding 
model at V( = 0.11 eV. The inset shows the V/ dependence of 
the top of the e' g band, A. 



Here, T c is decided under the condition that the eigen- 
value of the Eliashberg equation, A, is equal to one: 



k.n'.m 



V m G(k, ie n ,)<P(ie n ,) 



-is, 



,)V m D{ie n -ie n ,), (8) 



xG{- 

where we assumed a k-independent gap function 
4>{ie n ). Figure 4 shows the V{ dependence of T c . At 
(a SH , a BR ) = (0.2,0.2), T c increases as V{ decreases, 
whose change is much faster than that of the mass- 
renormalization factor. The obtained T c is high when 
the small pockets of the e' g band are present (V/ < 0.06 
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cV), reflecting a huge DOS at the Fermi level. It is note- 
worthy that T c keeps high values for a while even af- 
ter the small pockets disappear (V( > 0.06 cV), irre- 
spective of the fact that the DOS suddenly decreases 
by one-third at V/ = 0.06 eV. To find out the reason 
for this, we study the cases of (a SH , a BR ) = (0,0.2) 
and (a SH , ot BR ) = (0.2,0) for comparison. In the for- 
mer case, T c is low and independent of V/. In the lat- 
ter case, T c drops quickly when the small pockets dis- 
appear because shear phonons do not produce the at- 
tractive force between a\ g electrons. The reduction of 
T r for (a SH ,a BR ) = (0.2,0.2) is slower than that for 



(a SH , a BR ) = (0.2, 0), and T c is considerably higher than 
that for (a SH , a BR ) — (0, 0.2) even when the small pock- 
ets disappear. It is understood that the shear phonons 
support the SC and raise T c for a wide range of V{ . Be- 
cause the shear phonons enable cooper pairs to transit 
from the a\ g orbital to the e' g orbital, a sort of SK mech- 
anism works and T c increases. Indeed, the gap function 
4>e(uj) is finite only for the a l9 band for (a SH , a BR ) — 
(0,0.2) for V( > 0.06 eV; in contrast, 4>i{lo) is finite for 



both dig and e' bands when (a 



SH BR 



(0.2,0.2). By 



the analytical study for T c <C —A <C ujd, 17 we derive the 
expression T c m ujd cxp(— 1/A* ff ), where A* ff is given by 

2A^{ilog(^) + i} 



X* 



eff 



A* + 



l-AI{ilog(^) + i} 



(9) 



where X\ = 2 £z aigPaig , \* 2 - M 



A* ; 



26f 



-Z e i O e 



and AJ 



p e i . In the derivation of eq. 



(9), we put b 2 — and simplify the DOS as p aig (uj) = 
p aig , and p e > g (u) = p e ' g 8(A — u) per orbital. When 
|A| < ujd, the second term raises A* ff and T c . When 
\A\/lu d = 1/4 {V( ~ 0.09), \ log(^) + i ~ 1, and we es- 
timate that \\ = \* 2 k, 0.25, A^ w 0.50, and AJ w 0.13 in 
the present model. Thus, A* ff increases by 0.32 to a value 
of 0.57 owing to the second term, which is recognized as 
the valence-band SK effect. 
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Fig. 4. V/ dependence of critical temperature T c . Small hole 
pockets are absent when V( > 0.06. 



Finally, we comment on the possible origin of the 
anisotropy of the superconducting gap, which is sug- 



gested by the absence of the coherence peak in 1/T\T 
and by the specific heat measurement below T c . We have 
shown that the gap function of s-wave SC due to EPI 
became highly anisotropic when strong antifcrromag- 
netic (AF) fluctuation coexist, as realized in boron car- 
bide, (Y,Lu)Ni 2 B 2 C. 19 Thus, we expect that anisotropic 
s-wave state is realized in Nao.3sCo02 • J/H 2 because 
strong AF fluctuations are derived in the FLEX approx- 
imation. 12 This is an important future problem. We also 
comment on the effect of hydration on T c . The hydra- 
tion is expected to make V{ smaller since the thickness 
of C0O2 layer is reduced, which is actually suggested 
by NQR measurement. 20 Figure 4 tells that the reduc- 
tion of V/ due to hydration would explain why only hy- 
drated samples show SC. We stress that the e' g band in 
Nao.35Co02 • J/H2O does not intersect the Fermi level 
because the electronic heat coefficient 7 is almost un- 
changed by hydration. 1 Otherwise, 7 should increase to 
more than three times the experimental value by hydra- 
tion, according to the change in the DOS at the Fermi 
level. 12 In addition, the value of Knight shift is almost 
unchanged by hydration, which also indicate the absence 
of small hole pockets in hydrated samples. 5 Finally, we 
note that ref. 21 proposes a mechanism of s-wave SC due 
to the long-range Coulomb interactions. A mechanism of 
/-wave SC due to EPI and Coulomb interaction is also 
proposed. 22 

In summary, we discussed the possibility of s-wave 
SC for Na x Co02 based on a d-p model with a\ g and 
e' g bands, by considering EPI for breathing and shear 
phonons. The obtained T c is strongly enhanced when 
shear phonons are taken into account, in addition to the 
breathing phonons. The estimated EPI for both phonons 
are large enough to realize s-wave SC against the strong 
Coulomb interaction in Naa;Co02. We find that T c is still 
high even if small hole pockets are absent (V t > 0.6 eV), 
although the DOS at the Fermi level suddenly decreases 
to one-third of its value at V t > 0.6 eV. This enhance- 
ment of T c is brought about by the valence-band SK 
mechanism due to shear phonons. In the future, we will 
study the effect of Coulomb interaction to obtain a real- 
istic T c , and to explain the anisotropy of the supercon- 
ducting gap. 
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